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Abstract 

>: 

We discuss the double spin asymmetries in charmonium hadroproduction with nonzero 
transverse momenta at fixed target energies, yfs ~ 40 GeV, within the framework of the 
factorization approach. It is shown that the color octet contribution is dominant in the 
asymmetries. The size of the asymmetries and the projected statistical errors in a future 
option of HERA with longitudinally polarized protons scattering off a polarized target 
(HERA-iV) should allow to distinguish between different parametrizations for the polarized 
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gluon distribution in the proton. 
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The study of spin asymmetries in the production of heavy quarkonium states in polarized 



5_i ■ nucleon-nucleon collisions should provide important information about the proton spin structure. 



Heavy quark-antiquark production processes occur at small distances and the subprocess level 
cross sections can be calculated perturbatively. On one hand, charmonium production asym- 
metries are expected to be sensitive to the polarized gluon distribution function in the proton, 
since heavy quark systems are mainly produced in gluon-gluon fusion subprocesses. On the other 
hand, it is essential to investigate in more detail the heavy quark-antiquark pair hadronization 
phase. To this end, observation of charmonium production in polarized experiments is expected 
to provide additional tests for existing models. 

The double spin asymmetry in J/ ip production has been studied in the framework of the so- 
called color singlet model (CSM) |]J by Morii and collaborators [0. The CSM is a nonrelativistic 
model where the relative velocity, v, between the heavy constituents in a bound state is neglected. 
However, discrepancies between experimental data [3-10] and the CSM predictions indicate that 
0(v) corrections as well as other mechanisms of quarkonium production, which do not appear 
in the leading order in v, should be considered. Such an expansion of the quarkonium cross 
sections and decay widths has been realized in the last few years within the framework of the 



1 Talk given by A. Tkabladze at the Topical Workshop 'Deep Inelastic Scattering off Polarized Targets: Theory 
Meets Experiment', DESY-Zeuthen, September 1-5, 1997. 
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Factorization Approach (FA) based on Nonrelativistic QCD (NRQCD) ||11|| . Here the production 
cross section for a quarkonium state H in the process 



A + B -> H + X (1) 

can be written as 

a v = ^J dx i dx ifi/A{xi)fj/B{x2)o{ij -> H), (2) 

where fi/A is the distribution function of the parton i in the hadron A, and 

&(ij^lT)=ZC*[n](0\O H [n]\0), 

n 

the subprocess cross section, is separated into two parts: short distance coefficients, C y '[n], 
and long distance matrix elements, (0|C^[n] |0). The C lJ [n] is the production cross section of a 
heavy quark- ant iquark pair by fusion of partons i and j. It can be calculated in the framework 
of pQCD. The [n] state can be either a color singlet or a color octet state. The (0|O H [n]|0) 
describes the evolution of a quark-antiquark pair into a hadronic state. These matrix elements 
cannot be computed perturbatively, but the relative importance of long distance matrix elements 
in powers of velocity v can be estimated using the NRQCD velocity scaling rules ||12|| . Unlike 



the color singlet long distance matrix elements, each connected with the subsequent hadronic 
nonrelativistic wave function at the origin, color octet long distance matrix elements are unknown 
and should be extracted from experimental data. 

The new formalism implies that quark-antiquark color octet intermediate states are allowed to 
contribute to heavy quarkonium production and decay processes at higher order in the velocity 
expansion. Therefore, in the Color Octet Mechanism (COM) the complete structure of the 
quarkonium Fock space is taken into account while in the CSM only the dominant Fock state is 
considered, which consists of a color singlet quark-antiquark pair in a definite angular-momentum 
state of a final hadron (the leading term in the velocity expansion). 

The shape of the pr distribution of the 3 Si octet state production cross section indicates 
that J/ip and ip' production at large pt observed at the Tevatron (FNAL) can be explained in the 
FA JT0|, [13| . The color octet contribution to J/ip photoproduction has been analyzed in papers 
T3L IT5[. Recently, J/ip hadropro duct ion at fixed target energies has been studied by including 



the color octet mechanism |L6], [L7| [18[ . Large discrepancies between experimental data and the 
CSM predictions for the total cross section of J/ip hadroproduction were explained. The color 
octet contribution is dominant in J/ip hadroproduction at energies ^Js ~ 30 -=-60 GeV. The COM 
prediction for the ratio er(J '/ 'ip)dir I '&{J / V) — 0-6 is also in a good agreement with experimental 
data |I9J. 

Despite the obvious successes of the COM some problems remain unsolved. In particular, the 
theoretical predictions disagree with the J/ip and ip' polarization data at fixed target energies 



17] , |19|| and the COM prediction for the yield ratio of x& an d Xc2 states remains too low |L7 . 
These discrepancies indicate that higher twist corrections may give a significant contribution to 
low p T production of charmonium states and should be added to the color octet contributions 

20| . The color octet contribution underestimates the J ftp photoproduction cross section at large 
values of z (z = Ej/^/E y in the laboratory frame) [pT| . 



The NRQCD factorization approach implies universality, i.e. the values of long distance 
matrix elements extracted from the different experimental data must be the same. However, due 
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to the presently rather large theoretical uncertainties, the existing experimental data does not 
allow to check the FA universality and therefore to test the COM. 

This fact motivated us to look for other processes with less theoretical uncertainties to test 
the color octet mechanism. The observation of J/ if) asymmetries can be used for these purposes 
as well as measurements of the J /if) polarization in unpolarized hadron-hadron collisions and 
electroproduction |2(| |22|, |23| . 

In this talk we consider the possibility of extracting information about the gluon polarization 
in the nucleon through spin asymmetries in charmonium production. And conversely, we inves- 
tigate what can be learned from charmonium production asymmetries about color octet long 
distance matrix elements, if the polarized gluon distribution function should once be measured 
in other experiments. 

Finally, we will present the results of our calculation of the expected spin asymmetries in 
production of J /if} and Xcj states at HERA-iV, one of the future options of HERA 



an 



experiment utilizing an internal polarized nucleon target in the polarized HERA beam with 
energy 820 GeV would yield ^/s ~ 40 GeV. To avoid the uncertainties coming from higher twist 
subprocesses J20| we considered J/ if> production at large enough pt which can not be caused by 
internal motion of partons within the nucleon. 



2 Matrix Elements 

We consider the double spin asymmetry All in inclusive charmonium state production defined 

as 

H da(p+p + -> H) - do(p + p_ -> H) = EdAa/d 3 p 

LL[PP) da{p + p + — H) + da(p + p- -> H Edo/d?p ' 1 ' 

where p + (p-) stands for the sign of the helicity projection onto the proton momentum direction 
and H denotes the particular charmonium state. The production of each quarkonium state 
receives contributions from both color singlet and color octet quark-antiquark pair states. We 
consider here only the dominant sets of color octet states in the NRQCD velocity expansion for 
direct S and P state charmonium production. Hence we calculated the production asymmetries 
for the color octet states 3 -Po i 25 1 Sq and 3 Sf together with those for color singlet states. 

The number of color octet long distance matrix elements can be reduced using the NRQCD 
spin symmetry relations 

(O|0f( 8 P J )|O) = (2J+l)(0|Of( 3 P )|0), (4) 
(01^(^)10) = (2J+l)(0|Or( 3 Si)|0>, (5) 

which are accurate up to v 2 . 

After having utilized these relations we remain with the three independent color octet ma- 
trix elements (O^fSx)), (O^ ( 3 P )>, and (O^ C S Q )) which give the main contributions to 
the direct J /if) hadroproduction cross section. There is only one color octet matrix element 
(0g cl ( 3 Si)) which contributes to Xcj states production at lowest order. 

As was mentioned above, the values of the matrix elements extracted from the fixed target 
experiment's data contain large theoretical uncertainties. In our further calculations we will use 
the following values for the three main matrix elements extracted from J / if) production at CDF 
at the Tevatron using the GRV LO (1994) parametrization for unpolarized parton densities 
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(01^(^)10) = 
(0|O 8 J ^( 1 5 )|0) + ^(0jC% 



1.06 ±0.14+^ • lQ- 2 GeV 3 



J/4 (3 



mi 



Po)\0) 



3.9 ±1.15l};$- 10~ 2 GeV 3 



(6) 
(7) 



The second errors quoted in these expressions correspond to the variation of the factorization 
scale \i from 0.5 \jp\ ± ^ m c ^° ^\Jpt + 4m;?. From the errors indicated in (6) and (7) it is 
obvious that the variation of the renormalization and/or factorization scale also leads to large 
uncertainties when fitting the color octet parameters. 

The fit results for long distance color octet matrix elements can be affected also by higher 
v 2 corrections. In particular, large uncertainties emerge when the matrix element (Og^( 3 Si)) is 
extracted from the data on J / if) production at large pr at the Tevatron. In fitting the CDF data 
the energy of soft gluons emitted by a heavy quark-antiquark pair before transition into J j if) is 



usually neglected Jl3|. |22|. Ej| and, consequently, the fragmentation function of a gluon into J J if) 
on the scale 2m r has the form: 



D g ^(z,2m c 



7ra s (2m c 



Here the delta function implies that the J / if) carries the whole energy of the fragmenting gluon. 
If the nonzero energy of the soft gluons (of order m c v 2 ) is taken into account the fragmentation 
function becomes softer. Hence the realistic cross section is smaller at large and the fit value 
of the (Og^( 3 5'i)) matrix element appears about a factor of two larger [E7L ESI EM BO]. Therefore, 



care is required when using the value for this parameter extracted from CDF data at fixed target 
energies. Also, the uncertainty connected to the 'trigger bias' effect makes it impossible to use 
this value for testing the NRQCD universality. 

For indirect J /if) production via decay of the underlying if)' state and for production of Xcj 
states we use the following values fitted from CDF data: 



{OI (%)> + 



m 



(or i ( 3 s 1 )) 



9.8 • 10- 3 GeV 3 [13] 
0.46- 10- 2 Ge1/ 3 [|! 

1.8- 10~ 2 GeV 3 m 



(9) 
(10) 

(11) 



For the calculation of the expected asymmetries we assume that the first parameter is dominating 
in eq. (11). In this case the discrepancy between the COM predictions and experimental data 
on if)' polarization appears to be smallest [19 . 



3 Results and Discussion 

Asymmetries in production of quark-antiquark states. The characteristic value of the 
partonic x in the production of (cc) pairs can be obtained from the relation X\X 2 — (4m 2 + 
Pt)/S (~ 0.01 at HERA- A). This means that the typical values of x g i uon which can be 
probed by measuring the spin asymmetry in charmonium production is about x g i uon ~ 0.1. 
We used three parametrizations for polarized parton distribution functions (PDF) that are 
different in the region x ~ 0.1 to show the dependence of spin asymmetries in the produc- 
tion of various quark-antiquark pair states on the gluon polarization in the nucleon: the old 
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version of the Gehrmann-Stirling (GS) parametrization (set A) [31] (as example of a large 



gluon polarization) and the new version of the GS parametrizations in NLO and LO (both 
set A) |52| (as examples for moderate gluon polarizations peaking at different values of x g i uon ). 

In Fig.l the polarized gluon densities from these 
parametrizations are shown at Q 2 = 4 GeV (note 
that the polarized gluon distribution functions of 
NLO set B and LO set A practically coincide at this 
value of Q 2 ). As can be seen from Fig.l, the three 
chosen sets exhibit different values for the polarized 
gluon distribution function for the partonic x near 
0.1. We note that, although the calculations of sub- 
process level cross sections are performed in leading 
order, the NLO set of the parametrization was used 
to probe different shapes of the polarized gluon dis- 
tribution function. 

Fig. 2 shows the expected asymmetries for dif- 
ferent states of a heavy quark- ant iquark pair at 
y/s = 40 GeV (HERA- TV) . The asymmetries for 
and 3 S , j 8 ^ octet states are represented by solid 
and dotted lines, respectively. The dashed lines 
correspond to combined asymmetries of 3 Pj octet 
states, Sj=o 12 (2 J + l) 3 -^/ 8 ''- The dash-dotted lines 
show asymmetries of J/ip production in the CSM 
including J/ip production through decays of higher 
charmonium states {ip' and Xcj)- For unpolarized parton distribution functions we used the GRV 
LO parametrization 




Figure 1: Different possible polarized gluon 
distributions in the nucleon, used in this pa- 
per at Q 2 = 4 GeV: new GS parametrization 
NLO set A (solid line), LO set A (dashed 
line) and old GS parametrization (set A) 
(dash-dotted) El. 
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Figure 2: The expected asymmetries at HERA-V (y/s = 40 GeV) for the production of different color 
octet states. Solid lines represent asymmetries for the Sq state, dashed lines for the combination 
of 3 Pj 8 ^ states (see text), dash-dotted lines for the 3 s[ 8 ^ octet state, and dotted lines correspond to 
J/tp production in the CSM. The three figures base upon different input distributions: a) old GS 
parametrization (set A) new version of GS parametrization b) NLO set A, and c) LO set A [p^| . 



In all three parametrizations for the polarized PDF the gluon-gluon fusion gives the dominant 
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contribution to Act; quark-gluon subprocesses contribute only to about 10% and the contribution 
of quark- ant iquark annihilation subprocesses is less than 1%. 

As can be seen from Fig. 2, the expected asymmetries for all states strongly depend on the 
size of the polarized gluon distribution function in the region x g i uon ~ 0.1. 

Measurement of the inclusive J/ip asymmetry. In the inclusive case the kinematics 
of the 2 — > 2 subprocess cannot be reconstructed completely. Hence only indirect information 
on the gluon polarization in the nucleon can be obtained by measuring the spin asymmetry in 
inclusive J/tp production. 

Figure 3 shows the expected double spin asymmetries at HERA-iV energy as a function of J / ip 
transverse momentum. From now on we use the new GS parametrization, sets A,B, and C [p2| . 
For the mass of the charm quark the value m c = 1.48 GeV was taken and the parton distribution 
functions are evaluated on the factorization scale [l = \jpT + Am? c . The strong coupling constant 
is calculated by the one-loop formula with 4 active flavors {Kqcd = 200 MeV). 

Fig. 3a corresponds to the case when the first parameter in the combination (7), (Og^( l So)), 
is dominating, i.e. (0$^( 3 Po))=0. Fig. 3b corresponds to the opposite case, (Og // ^( 1 S'o))=0. In 
both cases we use for the third main parameter the value (0$^( 3 Si)} = 20 • 10~ 3 GeV 3 (the 
sensitivity of asymmetries to the value of this parameter will be considered below). Figures 3a 
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Figure 3: The expected asymmetries versus transverse momentum at \fs = 40 GeV for the NLO set 
A (solid lines), set B (dashed lines) and set C (dash-dotted lines) of the new GS parametrization p2| ; 
a) (O^ePo)) = 0, b) (O^CSo)) = 0. 

and 3b also display the expected statistical errors. The statistical error 5 All at HERA-iV can 
be estimated from |24| 

5A LL = 0.17/ ^pb)- (12) 

100% efficiency is assumed. This relation has been determined by assuming an integrated lumi- 
nosity of 240 pb~ l and beam and target polarizations Pb = 0.6 and P? = 0.8, respectively p4[ . 
The error bars are obtained by using integrated cross sections over bins Apr = 0.5 GeV (for the 
first three points) and Apx = 1 GeV (for the other two ones). The J/ip decay branching ratio 
into the e + e~ mode is also included. As can be seen from Figs. 3, the magnitude of asymmetries 
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and expected errors allows one to distinguish between different parametrizations of polarized 
parton distribution functions. 

Production asymmetry of Xci an d Xc2 states. In Figs.4a,b we show the production 
asymmetries of Xi an d Xi states in conjunction with the projected HERA-iV errors for the three 
different sets used for polarized PDF. Unlike direct J/ip production the production asymmetries 
of x states depend only on one color octet matrix element, namely ((9g cl ( 3 Si)). The value of 
this parameter was extracted from Tevatron data at large px values, i.e. from the region where 
the NRQCD factorization mechanism is valid, and higher twist effects are expected to be small, 
i.e. of the order of A/pr or even (A/pt) 2 - Hence the production asymmetries of x states do 
not contain additional free parameters and enable us to probe the polarized gluon distribution 
function with less theoretical uncertainties. We note that the branching ratios of xj decay into 
J/ijj plus photon are taken into account in the calculations of the expected errors. 
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Figure 4: The XcJ production asymmetries versus transverse momentum at \fs = 40 GeV for the NLO 
set A (solid lines), set B (dashed lines) and set C (dash-dotted lines) of the new GS parametrization 
H; a) Xd, b) X c2- 



J/ip plus jet production. Direct access to the ratio AG(x)/G(x) is possible only if the 
'other jet' is detected (photon+jet or J /ip+jet) f33fl , because the complete kinematics of the 2^2 
subprocess can be reconstructed if the away-side jet in the production of J/ip is measured, as 
well. If the values of the long distance matrix elements are established from experiments where 
the uncertainties connected with higher twist and other corrections are expected to be small (or 
negligible), the J/ip+jet production asymmetry at HERA-iV will serve as a good tool for the 
extraction of the AG(x)/G(x) value at x ~ 0.1 [|33| . 



Measurement of color octet matrix elements. As soon as the polarized gluon distri- 
bution function is extracted from other channels (e.g. photon+jet) at HERA-iV |33|] or RHIC 
31 1, a measurement of the J/ip asymmetry can be used for checking the NRQCD factorization 



scheme. 

For further calculations and analysis we choose the new GS parametrization for polarized 
PDF (NLO, set A) [j3^] . In figures 5a and 5b we present the expected asymmetries in J/ip 
production for two radical choices of (Og^ ( x So)) and (Cg^( 3 P ))- Figure 5a corresponds to the 
case when the parameter (0$^( 3 Pq)) tends to zero, i.e. (Og^^So)} is the dominating term in 
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Figure 5: The expected asymmetries at HERA- AT energy for different values of the long-distance 
parameter (0g^( 3 Si)}. The solid lines correspond to the value 20 • 10~ 3 GeV 3 , the dashed lines to 
10-KT 3 GeV 3 , and the dotted lines to 30-10~ 3 GeV 3 ; a) asymmetries for the case when (O^ ( 3 P )> = 0, 
b)for (O J 8 /lp ( 1 S ))=0. 



the combination (7). Figure 5b represents the other radical choice, when the second parameter 
(0g / V( 1 S'o)) j g zeTO \ s already mentioned above, due to the 'trigger bias' effect the value of the 
third main parameter, (0 8 ^ ( 3 Si)} , is underestimated about twice. To check the sensitivity of 

the expected asymmetries to the uncertainty caused 
by this effect, we used three different values for 
(Os 7 ^ 3 ^)), 10- lCT 3 GeV 3 from (6) (dashed lines in 
Figs.5), 20-lCr 3 GeV 3 (solid lines) and 30-lCr 3 GeV 3 
(dash-dotted lines). From comparison of Figs. 5a 
and b it can be seen that the asymmetries depend 
strongly on the relative values of the matrix elements 
(O^CSo)) and (O J 8 ^( 3 P )} and practically do not 
depend on the value of (0 8 ^( 3 Si)}. Hence, indepen- 
dently of the value of the latter parameter and, con- 
sequently, independently of the 'trigger bias' effect, a 
measurement of the J/ijj production asymmetry can 
10 " 1 1 10 be used for an extraction of the ratio of the first two 

parameters. 

Figure 6: The integrated double spin asym- In Fi S' 6 we P lot the expected integrated asym- 
metry versus the ratio of color octet parame- metr y (l- 5 < Pt < 3 GeV) versus the ratio a - 
ters, a = (O^ '(^K/ '(O^ \ 3 P )). Solid m^O^ C S Q )) / (O^ ( 3 P )) at HERA-iV energy for 
line- (0 J ^( 3 Si)} = 20 • 10~ 3 GeV 3 dashed the above presented three different values of the third 
line: 10 • 10~ 3 GeV 3 dotted line: 30 • 10~ 3 parameter. The expected statistical error for the in- 
GeV 3 . tegrated asymmetry is about 0.002, i.e. it will be 

possible to determine the ratio a within about ±15% 
if a is larger than unity. With a being smaller the sensitivity is becoming worse. 

Of course, this method requires that the polarized gluon distribution function is known 
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already and large enough to generate observable sizes of J/ip production asymmetries. 

One of the main parameters of the factorization approach is the mass of the charm quark. 
As was shown in [35[], the expected asymmetries are practically insensitive to the value of the 
charm quark mass. Therefore, the double spin asymmetry in J/ip production, unlike the cross 
section, should be free from uncertainties caused by the unknown mass of intermediate color 
octet states. We also note that the asymmetries do not depend strongly on the renormalization 
scale unlike the J / ip production cross section [^2] . 



Finally, we note that recently a possibility was shown to extract the ratio a of color octet 
long distance parameters from polarized J / ip electroproduction data with reasonable errors p3| . 



4 Conclusions 

We investigated the expected double spin asymmetries in heavy quarkonium hadroproduction 
in polarized proton proton collisions. To reduce the contribution from possible higher twist 
corrections |2(J we considered J/ip meson production at nonzero transverse momenta, px > 1.5 
GeV. Unlike the calculations of |36fl , where only the lowest order subprocesses, 2 — > 1, were taken 
into account, we considered J/ip production in 2 — > 2 subprocesses at large enough p T to avoid 
uncertainties coming from higher twist corrections |2(| . 

The size of the expected asymmetries in conjunction with the statistical errors at HERA-iV 
allow to distinguish between different parametrizations for polarized PDF (GS, set A, B and 
C p2|| ). If, on the other hand, those were already known measuring the asymmetry would 
give a possibility to extract information about the color octet long distance matrix elements 
(Og^( 1 S'o)) and (0^( 3 Po)), separately. This would deliver useful information to check the 
universality of the factorization scheme. 

A.T. acknowledges the partly support of this work by the Alexander von Humboldt Founda- 
tion. O.T. and A.T. wish to thank the Organizers of the Workshop and all the participants for 
the warmful atmosphere felt during their stay in Zeuthen. 
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